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Abstract 



o 

Q \ Bose-Einstein symmetrization can lead to correlations between out going 

\^ ' identical particles which reflect the space-time extent of the collision process. 

At LEP and LEPII these correlations have been studied as a function of the 



single variable Q = y— (Pi — P2) 2 - Assuming a simple form for the correlation 



C function the experiments find source radii dependent on the hadron mass. 



In this note, I point out that such effects can arise from purely kinematic 
£NJ ' considerations, although these are unlikely to explain the observed effects 

completely. 



I. INTRODUCTION 



Bose-Einstein correlations are of interest to the LEP community as a source of systematic 
error in W-mass determinations (for example |l|-f|]), and also to the Heavy Ion community 
(see |§,§]) as a valuable measurement tool for the size of the fireball created in heavy ion 
^ ! collisions (eg. [[?[ and refs therein). 

Due to limited statistics, measurements of BE enhancement at LEP are often made 

; 

in terms of a single variable Q = J — (pi — p 2 ) 2 where pi and p 2 are the 4-momenta of a 
pair of out going bosons. The enhancement is described by a correlation function which 
introduces an extra parameter, the HBT (Hanbury-Brown-Twiss [|J) radius. It has been 
observed II that the radii measured depend on the mass of the boson being used. In general 
the lighter bosons tend to give larger radii. In this section I describe the construction of 
the correlation function as a function of Q and its relation to the single particle Wigner 
distribution. In section 2 I use a simple Gaussian model for the Wigner density, and discuss 
the correlation function in various limits. In section 3 numerical results are presented for 
particles of various different masses and sources of different sizes. I conclude with a discussion 
of how far these results can go to explain the observed mass dependence of experimental 
Bose-Einstein parameters. 

The correlation function can be defined as 

c(Q) = ^ v (i.i) 
Pi(Q) 



l 



where p 2 (Q) is the density of particle pairs with invariant momentum separation Q, and 
Pi(Q) is the same density measured with respect to a reference distribution. The densities 
Pi and P2 can be defined in terms of the single- and two-particle momentum densities (Pi 
and P 2 respectively) according to 

Pi (q) = J d 3 Pl d 3 P2 Pi(pi)Pi(p 2 )«y(Q(pi,p 2 ) - Q) (1.2) 



P2(Q) = / ^ 3 Pi d 3 p 2 P 2 (p 1 ,p 2 )5(Q(p 1 ,p 2 ) - Q), 



;i.3) 



where the function Q(pi,p 2 ) is just J—(pi — £> 2 ) 2 - 

The single-particle momentum distribution is given in terms of the one-particle Wigner 
distribution S(x,p), 



^i(Pi) = J d A x S(x,p)\ p 2 =m 2. (1.4) 
It has long been known that under rather general assumptions the two-particle distribution 



can also be related to the one-particle Wigner density [JTO 
P 2 (pi,p 2 ) = Pi( P i)Pi(p 2 ) + 



d x S(x, K)e 



iq-x 



;i.5) 



where K — \{p\ + p 2 ) and q = p\ — p 2 . It is thus possible, given the source distribution, 
to compute the densities Pi and P 2 and hence pi and p 2 , then to construct explicitly the 
correlation function C(Q). The correlation function C(Q) may depend on the particle masses 
through the integration limits in equations (1.2) and (1.3), even when there is no explicit 
mass dependence in the source function S. In the next section I examine this implicit 
dependence for a Gaussian source. 



II. SIMPLE MODEL 

To illustrate the kinematical mass effects it is best to take a simple model. Consider an 
instantaneous, uncorrelated source, Gaussian in both space and momentum. The normalised 
source distribution is given by 

s( *' p,= (dp exp Hf) x 

The parameters Po an d Po are the typical size of the source in position and momentum 
space, respectively. The single and two particle distributions can now be obtained using the 
formulae of the previous section: 



2 



Pl(p) 

^2(P1,P 2 ) 



1 



exp 



|PT 

2Pn 2 



(2tt^)§ 
Pi(pi)Pi(p 2 ) + 

1 / IKI 2 



exp 



4P 2 



(2.2) 



(2.3) 



where K = |(pi + p 2 ) and q = (pi — P2). Note that in this simple example the fully 
differential particle distributions do not depend on the particle masses. 

Inserting Eqs. ( ^2| , ^3| ) into Eqs. ( p~72| , p~73|) and then constructing Eq. (|1.1| ) gives the 
expression for the correlation function 



C(Q) = 1 + 



/ d 3 pi d 3 p 2 exp 



|Pl+P2 

4P 2 



-P 2 | Pl - P2 | 2 ) 5(Q( P i,p 2 )-Q) 



fd 3 Pl d 3 p 2 exp (-^) exp (-Jg£) £(Q(pi,p 2 ) ~ Q) 



(2.4) 



Due to the spherical symmetry of the source most of the angular integrations are trivial. 
The (^-function constraint can be satisfied with the angle between Pl and p 2 ; define the 
cosine of this angle to be x. The integrations may then be re-written, 

d 3 Pl d 3 p 2 S(Q(p 1 ,p 2 ) -Q) — > p! dp x p 2 dp 2 Q(l - x)Q(l + x) (2.5) 

where Pi = \pi\, % = 1, 2. Constants independent of p\ and p 2 have been dropped as these 
cancel in the ratio ( |2.4| ). The variable function of p\ and p 2 is given by 



x 



(pi + m 2 )(p 2 + m 2 ) — m 2 
P1P2 



Q 2 



2pip 2 



(2.6) 



where Q(z) is the Heaviside (theta) function. Notice that mass dependence has now crept 
into C(Q) via the integration limits. 

The theta-f unctions in Eq. (|2.5|) give rise to a complicated integration region in the p\ —p 2 
plane which depends on the mass of the bosons. The integration region can be summarized 

as 



P2 > 



Pi 1 + 



Vi < Pi 1 + 



2m 2 



0^ 

2m 2 



Q 



\ 



1 



Pi_ 

m 2 



Q 2 ' 

Am 2 



+ Q 



\ 



p? N 

m 2 , 



1 + 



q 2 ' 

4m 2 



This takes on simple forms in the following limits 



Q 



m 



00 



P2 > 



4pi 



(the massless limit) 



Q 
m 

Q 
m 



Pi 1 < p 2 < Pi H 

V my \ m 

where pi,p 2 3> m. 

\pi -Q\<P2<Pi + Q 

where Pi,p 2 Cm 



(2.7) 
(2.8) 

(2.9) 

(2.10) 
(2.11) 
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Now the BE enhanced region of the p± — P2 plane lies on a band of thickness ~ Rq 1 
centered around p\ ~ P2, while the reference density lies largely in the region pi,p2 < O(Pq). 
Figures 1 and 2 show the enhanced region and the integration bounds for Q < m and Q > m. 
It is then clear that for any Q ^> m the enhanced region occupies a fraction 0Q.Ro.Po] 1 ) of 
the populated region. This factor must always be less than one by the uncertainty principle 
and so the strength of the correlation is reduced. In fact it is possible to show that in the 
massless limit, the intercept of the correlation function at Q = occurs at 



C(Q « 0) = 1 




e-a-atan- 1 ^)' 



e 3 j 

as f — > oo, (2.12) 



where £ = 2R P . 

The opposite limit Q -C m is qualitatively different. One can see from Eq. ([2.10|) that 



as Q — > with Q <C m, the integration region is squeezed into a narrow wedge along the 
diagonal of the p\ —p2 plane (the integration region in figure 1 is approaching this limit). If 
the inequality 

^ Q<^- (2.13) 
m Ro R Po 

is satisfied then the integration region only contains the enhanced region and the correlation 
function can approach the expected value of C(Q ~ 0) = 2. In fact in the limit that the 
mass is large compared to all other scales, the correlation function approaches the 'naive' 
form 

C(Q) = l + ex V (-Q 2 R 2 dass ) (2.14) 

where 

= Hi (l - pj^yi) ■ (2.15) 



III. NUMERICAL RESULTS 

In the case of an arbitrary mass the integrations must be done numerically. The corre- 
lation function shows different behaviour depending on the values of the mass and source 
parameters. The results of numerical evaluation of eqn.( |2.4|) are illustrated in figures (3,4). 
The typical momentum is set to the value 

P = 1 GeV, (3.1) 

in order to be comparable with the observed pion spectrum at the Z peak. This implies 
Rq > 0.1 fm which is consistent with all current estimates of the source size from BE 
Correlations. 
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Figure 3 shows the computed correlation function for a source of size 0.2 fm, and for 
three values of the mass roughly corresponding to the pion, kaon and A baryonQ masses. 
This corresponds to near the limit of the uncertainty principle (RoPo = \)- As implied in 
eq ( |2.12|) the correlation is more strongly suppressed in the region Q > m if RqPq is large. 
This is, in fact, quite general and the opposite is true; when -Ro-Po is close to the uncertainty 
limit the correlation function is only weakly suppressed. This is because when R$ 1 ~ Pq, 
the enhanced region extends over all the populated phase space (see discussion and figures 
from previous section) and averaging over the phase space causes no dilution. 

The opposite limit (RqPo large) is explored in figure 4. Here the source radius is R = 0.8 
fm which takes RqPo sufficiently far from the Heisenberg limit to give strong suppression. 
The effect is particularly noticeable for the pion correlation which becomes suppressed by 
a roughly a factor of 4 ~ RqPq in the region Q > m. This leads to typically non-Gaussian 
shapes for the pion (and other low mass) correlation functions. The suppression differentiates 
between the kaon and A correlations leading to a hierarchy of effective radii, R n > Rk > R\- 

From figures (3,4) one can see that smaller masses give characteristically sharper corre- 
lation functions as compared to heavier bosons. One can define an effective radius for the 
source (somewhat arbitrarily) in terms of the width at e _1 of the maximum of the correlation 
function. This results in the estimates shown in table 1 ( see [IT] and |I2j for pion and kaon 
measurements respectively). 

The reader may wonder if perhaps some intermediate source radius would provide a fit to 
the experimental radii. This is not the case. If one insists on fitting the source size measured 
by A correlations one is forced to a source of no larger than ~ 0.2 fm. The constraint that 
Po ~ 1.0 GeV now leaves -Ro-Po too small to provide the suppression necessary to fit the 
kaon and pion HBT radii. On the other hand, a source with RqPo > 1 leads to natural 
radii for pion and kaon correlations of order R T ~ m" 1 and Rk ~ m K which are larger 
than experimentally observed. An intermediate choice for Ro leads to distorted correlation 
functions which are partially suppressed for Q > m. These are difficult to parameterise in 
terms of a single radius. 

It is interesting to note that the correlation function for RqPq ^> 1 and m < P can be 
rather well approximated by the expression: 

C(Q) * 1 + e*p(- M { y±S^0^M \ (i±I±4 I , (3,) 



where 



y = • (3-3) 



in 



The suppression term multiplying the Gaussian tends to unity as Q — > 0, and approaches 
the value as Q — ► oo, where £ = 2RqPq as given in eqn. (|2.12|) . 



1 Here I plot positive correlations for the antisymmetric AA spin state for ease of comparision. 
The A is, of course, a fermion and after spin averaging displays anticorrelations, which lead to a 
suppression at small relative Q. This does not affect discussion of the HBT radius 
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IV. CONCLUSIONS 



I have shown that the averaging over phase-space to produce a correlation function 
of a single variable (Q) induces non-trivial, mass dependent distortions on the correlation 
function. It has been demonstrated that these distortions are such that correlations between 
smaller mass bosons are inherently shorter range in Q and lead to larger apparent radii. The 
effect is greatest away from the Heisenberg limit. 

For the physically reasonable value of P — 1.0 GeV, it is not possible to find a source 
radius R such that the kinematic averaging simultaneously fits the experimentally measured 
HBT radii in pion, kaon and A baryon correlations. As can be seen from table 1, away from 
RqPo ~ 1, averaging over phase-space can induce a mass dependence of the HBT radii of 
the same order of magnitude as experimentally observed (but not with the same absolute 
values) . 

Other effects are therefore necessary to bring theory into line with experimental data, 
but the effects of phase-space averaging should not be neglected. 
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TABLES 

TABLE I. Comparison of the radius parameters extracted from the simple Gaussian model (for 
Pq = 1 GeV, Rq = 0.2 fm and Rq = 0.8 fm) and experimentally measured values. 



Mass / GeV 


0.14 


0.5 


1.1 


(Source with Rq = 0.2 fm) Radius / fm 


0.32 


0.29 


0.22 


(source with Rq = 0.8 fm) Radius / fm 


2.8 


1.2 


0.9 


Experimental Radius / fm 


0.9 


0.6 


0.2 
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FIGURES 




j 

p, / GeV 



FIG. 1. A plot showing the allowed integration region given the parameters m = 0.14 GeV, 
Q = 0.05 GeV, Rq = 0.8 fm and Po = 1-0 GeV. The scatter shows schematically the enhanced part 
of the phasespace. For Q < m the integration is mainly over the enhanced region 




p, / GeV 

FIG. 2. A similar plot to fig 1, but with paramter Q = 0.3 GeV. With Q > m the allowed 
integration region is not restricted to the enhanced region and the correlation is diluted. 
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FIG. 3. Correlation functions computed for boson masses ~ 0.14 GeV, mx ~ 0.5 GeV 
and ttt-a ~ 1.1 GeV. The dashed curve shows the correlation function in the extreme heavy mass 
limit which corresponds to the naive expectation. The correlation functions are computed from a 
Gaussian source of radius Rq = 0.2 fm and typical momentum Pq = 1.0 GeV, Q = \/—(pi — P2) 2 
- the invariant momentum separation. 
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FIG. 4. Correlation functions computed for boson masses ~ 0.14 GeV, mx ~ 0.5 GeV 
and m\ ~ 1.1 GeV. The dashed curve shows the correlation function in the extreme heavy mass 
limit which corresponds to the naive expectation. The correlation functions are computed from a 
Gaussian source of radius i?o = 0.8 fm and typical momentum Pq = 1.0 GeV, Q = y/—(pi — P2) 2 
- the invariant momentum separation. 
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